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EFFECT  OF  THE  ACCELERATION  OF  ELONGATED  BODIES  OF  REVOLUTION 
UPON  THE  RESISTANCE  LN  A  COMPRESSIBLE  GAS* 
By  F.  I.  Frankl 


The  problem  of  the  motion  of  an  elongated  "body  of  revolution  in 
an  Incompressible  fluid  may,  as  is  known,  "be  solved  approximately  with 
the  aid  of  the  distribution  of  sources  along  the  axis  of  the  body.   In 
determining  the  velocity  field,  the  question  of  whether  the  body  moves 
uniformly  or  with  an  acceleration  is  no  factor  in  the  problem.  The 
presence  of  acceleration  must  be  taken  into  account  in  determining  the 
pressures  acting  on  the  body.  The  resistance  of  the  body  arising  from 
the  accelerated  motion  may  be  computed  either  directly  on  the  basis  of 
these  pressures  or  with  the  aid  of  the  so— called  associated  masses 
(inertia  coefficients).  A  different  condition  holds  in  the  case  of  the 
motion  of  bodies  in  a  compressible  gas.   In  this  case  the  finite 
velocity  of  sound  must  be  taken  into  account.   If  it  is  assumed  that 
the  body  produces  in  the  flow  only  small  disturbances,  the  velocity 
potential  cp  satisfies  the  wave  equation: 

a2^p  =  &  (1) 

where  a  is  the  velocity  of  sound.  The  method  of  sources  still  remains 
applicable  but  in  computing  the  effect  of  the  sources  it  is  necessary 
to  make  use  of  the  retarded  potential. 

We  introduce  two  systems  of  coordinates,  one  a  fixed  system  xr, 
the  other  moving  with  the  body  £r,  where  x  and  £   are,  respectively, 
the  coordinates  along  the  axis  and  r  is  the  distance  from  the  axis. 
The  coordinate   I   is  computed  from  the  nose  of  the  body. 

Let  r  =  r( !)   be  the  equation  of  the  body  of  revolution  (0  <  I  <  l) . 
The  position  of  the  body  is  characterized  by  the  inequalities 
— f(t)  <  x  <  Z  — f(t),  where  f(t)   is  a  given  function  of  the  time 
characterizing  the  motion  of  the  body.  The  coordinates  x  and  I  are 
connected  by  the  equation  x  =  |  —  f  (t)  . 

*"0  Vliianii  uskoreniya  na  soprotivlenie  pri  dvizhenii  prodolgovatykh 
tel  vrashcheniya  v  gazakh."  Prikladnaya  Matematika  i  Mekhanika,  Vol.  10, 
no.  k,   191*6,  pp.  521-521*. 
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Let  q.  =  q/x*  ,t' )   "be  the  value  of  the  linear  density  of  the  source 

at  the  point  x'   at  the  instant  of  time  t'   and  E  =  \J  (x  —  x' ) 2  +  r2 
he  the  distance  between  the  source  and  the  point  considered.  The 
potential  of  the  disturbance  velocities  is  then 


^^  --^jh^^-l)^ 


where  the  integral  is  taken  between  the  limits  determined  by  the 
inequalities 


_f  (t_D<x.<z_f  (t-f) 


In  order  to  establish  the  range  of  integration,  we  must  consider 
in  the  plane  x't'   a  strip  -f(t')  <  x'  <  I   -  f  (f )  between  the 
lines  of  motion  of  the  front  and  rear  points  of  the  body  (fig.  l)  and 
the  part  of  the  semihyperbola  in  this  strip. 


a 


-i/U'  -x)2^-K2 


The  abscissas  of  the  part  of  the  semihyperbola  considered  form  the 
range  of  integration  (fig.  l) .  As  may  be  seen  from  the  figure,  the 
range  of  integration  may  consist  of  one  or  several  segments.  The  ends 
of  the*  corresponding  arcs  of  the  hyperbola  may  lie  either  on  the  line 
of  motion  of  the  nose  or  on  the  line  of  motion  of  the  base.   In 
particular,  in  the  motion  of  the  body  with  supersonic  velocity  it  is 
possible  that  within  the  strip  considered  there  lies  an  entire  arc  of 
the  hyperbola,  both  ends  of  which  lie  on  the  line  of  motion  of  the 
nose  (as  will  be  the  case,  for  example,  for  a  uniform  supersonic  motion 
or  for  a  motion  approximating  this  type) . 

The  function  q(x' ,t')  must  be  determined  from  the  integral 
equation  expressing  the  fact  that  the  normal  component  of  the  velocity 
of  the  gas  on  the  surface  of  the  body  is  equal  to  the  normal  component 
of  the  motion  of  the  body: 


cki     kn   Sn. 


M^-l) 


cbc.  =-_^L=f«(t)     (2) 
r=r       \[l  +   [r«(ij]  d 


where  d/cta  denotes  differentiation  along  the  outer  normal  and 
%    =x  +  f(t).   In  what  follows  we  consider  the  derivative  br/bt,     as 
small  and  shall  give  an  approximate  solution  of  equation  (2).  For  the 
assumption  made,   oXp/dn  may  be  replaced  by  ckp/oT.  We  assume  for 
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simplicity  that  the  range  of  integration  reduces  to  the  segment   (x^Xg) 
Equation  (2)  then  assumes  the  form 


»2  V2     a/  or    Ri  V  1'    a/  ^r      v  '  d| 


«L 


dx« 


r=r 


(3) 


If  r/l  and  dr/d|  are  considered  small  quantities  of  the  first  order 
and  their  higher  degrees  or  derivatives  are  neglected,  an  approximate 
solution  of  equation  (3)  will  be1 


q(x,t)  =  2*f»(t)r(|)  g 


(*) 


We  shall  show  that  the  error  in  substituting  this  expression  in 
equation  (3)  is  a  magnitude  of  the  second  order  of  smallness.  For  this 
prupose  it  is  convenient  to  write  e r*  for  r  where  r*  =  r*(|)   is  a 
variable  quantity  of  the  order  of  unity  relative  to  I     and  6   is  a 
small  constant  number.  The  degree  of   e  in  each  expression  then  shows 
its  order  of  smallness.  We  consider  the  order  of  the  terms  of  the  left 
side  of  equation  (3).  The  terms  outside  the  integral  on  the  left  of 
equation  (3)  containing  q  as  a  factor  for  finite  terms  (in  the  general 
case)  will  be  of  the  second  order  of  smallness.   In  the  same  way  the 
expression 


dq  djc' 


E< 


<Me2r 


^-<M*e2 


E' 


(5) 


is  of  second— order  smallness.  We  write  the  remaining  term  in  the  form 


r 


i  c  &*  -  ?) w = &m  f  s  *  s  r  **■*  -  iu't,] 


dx» 
E3 


q(i«,t  -|)  -q(x«,t) 


dx* 
e3 


(6) 


-^his  solution  (k)  was  first  used  by  Karman  (reference  2)  in  the  case 
of  a  uniform  motion  (in  particular,  for  obtaining  the  approximate  value  of 
the  wave  resistance) . 
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As  a  result  of  the  computation  of  the  first  integral  of  this  expression 
we  obtain 


da* 

E3 


=r*(£) 


2  + 


0(e2)] 


(7) 


For  the  estimates  of  the  differences  in  the  brackets  in  the  second 
and  third  integrals  of  expression  (6)  we  have,  respectively 

a(i»,t)  -  q(x,t)  =  0(e2)E,     q(x',t)  -  q(xl,t  -  E/a)  =  0(e2)E     (8) 

Hence  the  estimates  of  integrals  are  of  the  form 

1x2 


0(e2)r 


'*! 


ax1 
E2 


<  «0(e^) 


(9) 


so  that  these  integrals  are  likewise  magnitudes  of  the  second  order  of 
smallness .  Thus 


kit   dr 


1X2  r 


J*l 


K*.*-S 


r=r 


dx«  =  f'(t)  ||  +  0(e2) 


(10) 


as  was  required  to  be  proven. 

The  expression  for  the  velocity  potential  after  substituting  in 
equation  (4)  becomes 


q>(x,t)  =  -  I  [|  ft  (t  -  I)  7(V)   dx»    (j(0   =  x   ff)        (11) 


For  computing  the  field  of  pressures  we  have  the  generalized  formula 
of  Bernoulli 


P  -  P     vd       dcp 

p  2  "  at 


(12) 


where  p  is  the  pressure  at  the  given  point,  p  the  pressure  in  the 
undisturbed  region,   p  the  density  and,  w  the  modulus  of  the  velocity. 
Neglecting  magnitudes  of  the  second  order  w2/2  and  taking  account  of 
the  equations 


b-|j 


ot- 


=  ft 


\   a/ 


(13) 
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we  obtain 
P  -  P 


1 
2~ 


'xl 


2R 


(t  -|)12  7*(V)    +  f"  (i 


-|)  7(1') 


dx1 
R 


•=X2 


S12 

~oT 


(HO 


In  the  particular  case  where  the  pressure  is  required  on  the  surface 
of  a  projectile,  the  motion  of  which  approximates  to  uniform  supersonic 
motion,  the  point  Xg  lies  on  the  line  of  motion  of  the  nose  of  the 

projectile  and  the  last  term  on  the  right  of  equation  (l^)  is  equal 
to  zero.  We  obtain 


P  -  P   1 
P     2 


,x2fr 


/xi 


f '  t 


-!) 


7*(V)   +  f " 


M) 


dx' 
R 


^      (15) 


where  x-,   and  Xg  are  determined  from  the  equation 


%2  =   f  (*  -kKxl,2-¥  +  **) 


In  the  case  of  uniform  motion  x  =  |'  —  vt  +  c  and  equations  (15) 
and  (l6)  become 


IX  * 


(16) 


p  -  P 
p 


(17) 


JX  * 

;il 


The  additional  pressure  produced  by  the  acceleration  is  therefore  given  by 
the  equation 


Sp 


,x2 


=  1    "f»(t-  |)7(5') 


dx' 
R 


>x. 


Jx]_ 


f,(*-|)2-  [r^)]2^') 


dx' 
R 


[f(t)]a  jj^V'fe')  -  r'(5«3^+  £**&<*>]  2  7-(iM 


dx' 
R 


.  [f'(t)]2  r'(i') 


2  v.  r  n  ^  2*1 

R 


(18) 
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where  £**   is  determined  from  the  equation 

xt  =  g*i  _  ft(t)t»  +  f'(t)t  -  f(t)  (19) 

From  the  structure  of  equation  (l8)  it  is  easily  seen  that  the 
relative  increase  of  the  pressure  arising  from  the  acceleration  has,  for 
velocities  comparable  vith  the  velocity  of  sound,  an  order  of  magni- 
tude b2/v2  where  h  =  f"(t)   is  the  acceleration.  Hence  for  rocket 
projectiles  of  the  usual  dimensions,  for  accelerations  not  higher 
than  1000  meters  per  second2  and  velocities  comparahle  with  the  velocity 
of  sound,  the  added  pressures  arising  from  the  acceleration  are 
negligibly  small. 

This  is  confirmed  also  by  the  following  illustrative  computation. 
The  length  of  the  war  head  portion  of  the  shell  is  equal  to  I   =  0.25  meter 
and  the  maximum  radius  r  =  0.07  meter.  The  generator  of  the  head  is  an 
arc  cf  a  parabola  which  smoothly  goes  over  into  the  cylindrical  part  so 
that  it  is  given  by  the  equation  r  =  O.56E,  —  1.12|2  (scale  in  meters). 
The  motion  is  one  of  uniform  acceleration  with  an  acceleration 
of  1000  meters  per  second2  so  that  the  line  of  motion  of  the 
projectile  x  =  -5000t2   (t  >  0),  x  =  0   (t  <  0) . 

The  pressure  distribution  over  the  war  head  is  found  at  the  instant 
of  time  t  =  0.5  second  such  that  the  velocity  is  equal  to  v  =  500  meters 
per  second.   The  air  density  p  =  0.125  kilograms~lper  second2.   There 
were  obtained  for  the  pressure  distributions  (with  account  taken  of  the 
acceleration)  and  the  additional  pressures  6p  produced  by  the 
acceleration  the  following  values: 

i   =  0.05  0.10  0.15  0.20     I    =  0.05  0.10  0.15  0.20   (meters) 
p  -  p  =  8550,71^8,5970,3588;  6p  =  -5.0  -5.5  -11.0  -16.0  (kilograms  per  meter2) 

The  wave  resistance  and  the  corresponding  added  resistance  were  then 
obtained  as 

Q  =  2rt  /   r  J|  (p  -  p)  d£  =  100,    SQ  =  2*  /   r  ^8p  d£  =  0.13   (kilograms) 
J°     '  i0  l  (20) 
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The  added  pressures  Sp  arising  from  the  acceleration  were  computed  from 
the  approximate  equation  (18)  which  for  the  example  given  has  the  form 
r 


Sp 


=  £b 


/    22l|Ildx'   -2^     /    2  7'(i»)   dx«    +(^)2     /    2  7"  (&»)   R  dx' 


b7«(0) 


tg  0x2       t]_  6x-j_ 


R2 


where 


5xt    = 


i  -  *i  - 


xx*  s 


2W 


bt  *1 
a 


El 


-  1 


6x2  =  X2  —  X2* 


0 


bt  x  -  x2 
a    '    E2 


-  1 


(21) 


(22) 


Translated  by  S.  Eeiss 
National  Advisory  Committee 
for  Aeronautics 
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